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Type I planetary nebulae (PNe) have high He/H and N/0 ratios and are 
thought to be descendants of stars with initial masses of ~ 3 - 8Mq. These 
characteristics indicate that the progenitor stars experienced proton-capture nu- 
cleosynthesis at the base of the convective envelope, in addition to the slow 
neutron capture process operating in the He-shell (the s-process). We compare 
the predicted abundances of elements up to Sr from models of intermediate-mass 
asymptotic giant branch (AGB) stars to measured abundances in Type I PNe. 
In particular, we compare predictions and observations for the light trans-iron 
elements Se and Kr, in order to constrain convective mixing and the s-process 
in these stars. A partial mixing zone is included in selected models to explore 
the effect of a ^^C pocket on the s-process yields. The solar-metallicity mod- 
els produce enrichments of [(Se, Kr)/Fe] < 0.6, consistent with Galactic Type 
I PNe where the observed enhancements are typically ^0.3 dex, while lower 
metallicity models predict larger enrichments of C, N, Se, and Kr. O destruction 
occurs in the most massive models but it is not efficient enough to account for 
the > 0.3 dex O depletions observed in some Type I PNe. It is not possible to 
reach firm conclusions regarding the neutron source operating in massive AGB 
stars from Se and Kr abundances in Type I PNe; abundances for more s-process 
elements may help to distinguish between the two neutron sources. We predict 
that only the most massive (M > SM©) models would evolve into Type I PNe, 
indicating that extra-mixing processes are active in lower-mass stars {3-4Mq), if 
these stars are to evolve into Type I PNe. 

Subject headings: stars: AGB and post- AGB stars — planetary nebulae: general 
— nuclear reactions, nucleosynthesis, abundances 



Introduction 



After the thermally-pulsing AGB (TP-AGB) phase is terminated, low to intermediate 
mass stars (~ 0.8 to 8Mq) evolve to become post-AGB stars and possibly planetary nebulae 
(PNe), if the ejected envelopes have sufficient time to become ionized by the hot central 
stars before dissipat ing into the in terst ellar mediuna . For recent reviews of TP-AGB and 
post-AGB stars see iHerwigl (120051 ) and Ivan Winckell (120031 ). respectively. The illuminated 
PN is comprised of material from the deep convective envelope, hence nebular abundances 
should reveal information about the efficiency of mixing events and chemical processing that 
took place d uring previous evolutionary phases, in addition to the initial composition of the 
parent star (IDopita et al.lll997l : iKarakas fc Lattanzidl2003al ). 



- 3 - 



Briefly, during the TP-AGB phase the He-burning shell becomes thermally unstable 
every 10^ years or so, depending on the core mass. The energy from the thermal pulse 
(TP) drives a convective zone in the He-rich intershell, that mixes the products of He- 
nucleosynthesis within this region. The energy provided by the TP expands the whole star, 
pushing the H-shell out to cooler regions where it is almost extinguished, and subsequently 
allowing the convective envelope to move inwards (in mass) to regions previously mixed 
by the flash-driven convective zone. This inward movement of the convective envelope is 
known as the third dredge-up (TDU), and is responsible for enriching the surface in ^^C 
and other products of He-burning, as well as heavy elements produced by the s process. 
Following the TDU, the star contracts and the H-shell is re-ignited, providing most of the 
surface luminosity for the next interpulse period. In intermediate-mass AGB stars with 
initial masses > 4:Mq, the base of the convective envelope can dip into the top of the H- 
burning s hell, causing proton-c apture nucleosynthesis to occur there (hot bottom burning, 
HBB; see Lattanzio et al. 1996. for a review). 



Abundances can be obtained from PN spe ctra for a number of elements including C, N, 



O, S, CI and the noble gases He, Ne, and Ar f|Aller fc Czvzak 



19941 : iDopita et all 119971 : IStanghelhni et al.ll2000l : iLeisv fc Dennefeld 



19831: 



Kingsburgh &: Barlow 



ments synthesized by t he s process, including Ge, Se, Kr, Xe, Ba, and Br (IPeauignqt fc Baluteau 



1994 lDinersteinll200ll : [Sterling et al.ll2002l : ISharpee et al.ll2007l : ISterhng et al.ll2007l : [Sterling fc Dinerstein 



1^ 

20061): and heavy ele- 



20081 . hereafter SD08). Whereas the metallicity of a star is usually taken from the abundance 
of iron, for ionized nebulae oxygen is generally the primary index of metallicity. However, 
O is not an ideal metallicity index for PNe, because it can be destroyed by CNO cycling 
and / or synthesized during hel ium burning. Consequently, Ar is sometimes used in place 
of O (ILeisy fc Dennefeldl l2006l . SD08), under the assumption that it remains unchanged by 
AGB nucleosynthesis. Zinc is a better indicator of iron-group element ab undances in nebulae 
than Fe itself. 



Sembach et al 



2001 



jecause Zn does not co ndense into dust as easily as Fe (lYork fc Jura 



Dinerstein 



9951: 



Welty et al.l 119991 ) and can be observed in PNe (IDinerstein fc Geballe 



1982 



20041 ). However, Zn is at the beginning of the s-process chain and might be 
produced to some extent in AGB stars. In order to interpret the elemental abundances in 
PNe it is important to determine which elements are produced, destroyed or left unaltered 
by AGB nucleosynthesis. 

The set of elements heavier than Fe that can be detected in PNe include all even- 
numbered elements from Z = 30 - 36, which lie at the light-element end of the s-process 
distribution, j u st be low the first s-process peak at Sr, Y, and Zr (Z = 38 - 40). In 
Karakas et al.l (20071) we at tempted to use Ge abundances derived from a small sample 
of PNe (ISterling et al.ll2002l ) to try and constrain the amount of TDU mixing after the final 
few TPs. There a hint was found that efficient TDU at the tip of the AGB may be required 
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in order to match the abundances; however, stronger conclusions were precluded by model 
uncertainties, the small PNe sample size, and the large uncertainties of the Ge abundance 
determinations. In the present study we proceed to analyze Se and Kr abundances from a 
much larger sample of 120 Galactic PNe from SD08, where the s-process enrichments were 
investigated for correlations with PN morphology and other nebular and stellar characteris- 
tics, as outlined in detail in ^ We also use PNe abundances of C, N, and O along with the 
abundances of post-AGB stars for comparison to the stellar models. 

The composition of Type I PNe provides important constraints on nucleosynthesis and 
mixing in intermediate-mass AGB stars. These PNe tend to exhibit bipolar morpholo- 



(Stanehellini et al. 


20 


06 




~3 - 8 M0 f 


Peimbert 


1978 



AGB stars, on the other hand, are defined as h aving initial mass es 



3M(7) a nd tend to 



show carbon and s-process element enrichments (IBusso et al.ll200ll : iTravaglio et al.ll2004i ). In 
the Milky Way, the spatia l distribution and kinema tics of Type I PNe indicate that they are 



a young population (e.g., ICorradi fc Schwardll995l ). This supports their identification with 



intermediate-mass stars, which have relatively short evolutionary lifetimes. In the Galaxy, 
intermediate-mass AGB stars are difficult to identify owing to a lack of reliable distances. 
This has resulted in a paucity of observational evidence for constraining stellar models, 
and especially the e fficiency of the TDU in intermediate-rnass AGB stars, which is still 



not w el l determined (IKarakas et al.l 120021 : 1 Ventura et al.l 120021 : lHerwig||2004j : IStancliffe et al. 
2OO4J ). iGarcia-Hernandez et al. J2006h identified several Galactic intermediate-mass AGB 
stars within a sample of OH/IR stars (i.e., bright 0-rich giants with large infrared excesses). 
The large enhancements of Rb found by these authors, combined with the fact that these 
stars are 0-rich, support the prediction that HBB and an efficient TDU has occurred in 
these stars. At lower metallicities, constraints are provided by observations of luminous O- 
rich AGB stars in the Large and Small Magellanic Clouds (LMC and SMC, respec t ively) that 



are r ich in Li and s-process elements ( IWood et al.lll983l ; ISmith &: Lambertlll989l ; iPlez et al. 
1993h . 



We begin in §2] with a summary of observational results for n-capture elements in PNe. 
We review the numerical method and present the stellar models in §31 and outline the results 
in §11 In §5l we discuss our findings and their implications, and in §61 we summarize our 
conclusions. 
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2. Observational Constraints 

The low cosmic abundances of n-capture elements cause their spectroscopic features 
to be very weak. For that reason, trans-iron elen ients were not identified in the spectrum 
of a PN until 1994 (jPequignot fc Baluteaul Il994l ). However, in the last few years there 
has been considerable progress in the detection of emission lines from and derivation of 



abundances for trans-iron elements in PNe (e.g., iDinerstein fc Geballdl200ll : [Sterling et al. 
200?! : ISharpee et allboOTl . SD08). 



The most comprehensive study is that of SD08, who conducted the first large-scale 
survey of n-capture elements in PNe. They determined Kr and Se abundances in 120 PNe 
from the n e ar-inf rared emission lines [Kr III] 2.199 fim and [Se IV] 2.287 /im identified by 
DinersteinI (120011 ). and corrected for the abundances of unobserved Se and Kr ions usin g 
analytical formulae derived from a grid of photo-ionization models (ISterling et al.l 120071 ). 
This survey increased the number of PNe with determined n-capture element abundances 
by nearly an order of magnitude, enabling a detailed study of s-process enrichments in PNe 
as a population, as well as an assessment of enrichment patterns in different classes of PNe. 
In order to normalize the abundances of the s-process products to the initial metallicity of 
each PN, the Se and Kr abundances were compared to those of O and Ar reported in the 
literature. Oxygen was utilized as a reference element for PNe with low-mass progenitors, 
since its abundance is generally the most reliably determined in PNe, while Ar was used 
for objects with higher-mass central stars that may have experienced O destruction during 
HBB. Overall, 41 of the 94 PNe with derived Se and/or Kr abundances or meaningful upper 
limits were found to be enriched in these elements, with the average [Se/(0, Ar)]= 0.31 and 
[Kr/(0, Ar)]= 0.9^1. SD08 interpreted these enrichments as evidence for in situ s-process 
nucleosynthesis and dredge-up in PN progenitor stars. 

SD08 observed significant numbers of PNe with intermediate-mass progenitor stars, as 
evidenced by their Type I compositions and/or bipolar morphologies (29 and 28 objects 
respectively; see their Table 15). They found that Type I and bipolar PNe exhibit smaller 
Se and Kr enrichments on average than other PNe. In Figured] we display [Se/(0, Ar)] and 
[Kr/(0, Ar)], separated into 0.1 dex bins, for non-Type I, Type I, and bipolar PNe. The 
distribution of Se and Kr enrichments is clearly skewed to smaller values for Type I PNe 
than other objects, and to a lesser extent the same is true of bipolar PNe. Notably, all of 
the Type I PNe observed by SD08 exhibit modest (less than a factor of two, or 0.3 dex) or 
no Se and Kr enrichments. These values are significantly lower than the average values for 



^We use the notation [Se/(0, Ar)] and [Kr/(0, Ar)] of SD08 to emphasize the use of different reference 
elements for PNe with different progenitor masses, where [X/Y] = logio(X/Y) — logio(X/Y)Q. 
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non-Type I PNe: 0.36 dex for [Se/(0, Ar)], and 1.02 dex for [Kr/(0, Ar)]. In contrast, while 
some bipolar PNe are among the objects with the lowest [Se/(0, Ar)] and [Kr/(0, Ar)], 
others display enrichments approaching a factor of 10. 

Figure [2] shows [Se/(0, Ar)] and [Kr/(0, Ar)] as a function of the logarithmic He/H, 
N/0, and C/0 abundances of PNe from the SD08 sample, including upper limits (open 
symbols) for objects without detected Se and/or Kr emission. This figure not only illustrates 
the lower average Se and Kr abundances of Type I PNe relative to other objects, but also that 
upper limits to the Se and Kr abundances in many Type I PNe allow for at most marginal 
s-process enrichments. SD08 found that s-process enrichments are positively correlated with 
C/ O in non-Type I PNe ( bottom pariels of Figure El), a result that agrees with studie s of AGB 
JSmith fc Lambertlll990l : kbia et al.l[2002[ ) and post- AGB stars Jvan Winckell boosh . Type I 
PNe often exhibit low C/0 ratios (< 1), in agreement with predictions from intermediate- 
mass stars with HBB. HBB converts dredged-up ^^C int o ^'^N, thereby preventing or delaying 
the formation of C-rich AGB stars (IFrost et al.lll998l ). On the other hand, some bipolar 
PNe display significant enrichments of both C and s-process nuclei, as well as non-Type I 
compositions. These objects may be descendants of low- mass binary star systems in whic h 
the bipolar morphology is a result of binary interactions (jSokerlll997l : iBalick fc Frankll2002l ). 
This suggests that some bipolar PNe are not descendants of intermediate-mass stars, and 
that morphology is not as reliable an indicator of progenitor mass as is chemical composition. 



Sharpee et al.l (120071 ) derived n-capture element abundances in five PNe from high reso- 



lution optical spectra. They identified emission lines of several n-capture elements, including 
Br, Kr, Xe, Rb, Ba, and Pb in each of these objects, and derived abundances for Br, Kr, 
and Xe. They found that Kr and Xe are enriched in three of the five observed PNe, by 0.3 
to 0.9 dex. The two objects in their sample which do not exhibit 5-process enrichments, 
NGC 244 and IC 2501, are both Type I PNe according to the classification scheme of 
PeimbertI (119781 ). The [Kr/Ar] and [Xe/Ar] values of these two objects are consistent within 
the uncertainties with the solar values, except for Xe in NGC 2440, which may be mildly 
subsolar. These results for [Kr/Ar] are consistent with the findings of SD08 for their larger 
sample of Type I PNe. 

The near-infrared [Kr III] and [Se IV] lines utilized by SD08 have also been detected 



in two PNe belonging to metal-poor Milky Way dwarf satellite galaxies (IWood et al.l 12006 



Dinerstein et al., in preparation). Such observations can provide important constraints 
on AGB models and nucleosynthetic yields for subsolar metallicities. One of the observed 
objects, LMC SMP 62, is N-rich (N/ O ~ 0.5) and hence considered a Type I PN in the context 
of the LMC, w here the threshold value of N/0 for Type I objects is takeri to be lower than for 
Galactic PNe JPopita fc Meatheringhamlll99ll : lleisv fc Dennefeldll2006l ). While the [Kr III] 
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line is not detected (and does not provide a useful abundance limit), the [Se IV] line yields a 
Se abundance ~ 15% solar (Dinerstein et al., in preparation). Not only is Se not enriched, 
it is actually deficient relative to the a-elements: the abundances of O, Ne, S, an d Ar in 
SMP 62 are - 30-40% solar Jleisv fc Dennefeldl bood : ISernard-Salas et all boOsl lR. This 
is reminiscent of the abundance pattern in LMC red giants, where light g-p rocess elements 



such as Y and Zr are marked ly deficient relative to iron and some a elements (jPompeia et al. 



20061 : iMucciarelli et al.ll2008l ). In contrast. Hen 2-436, a PN in the Sagittarius dSph galaxy. 



displays Se and Kr abundances as high as those for some of the most enriched Galactic 
objects shown in Figure [H However, Hen 2-436 is C-rich rather than N- rich, and therefore is 
a non-Type I PN, with an estimated initial stellar mass of ~ 1.2 Mq (jPudziak et al.ll2000l : 



Ziilstra et al.ll2006f ). 



Zinc is detectable in P Ne via the [Zn IV] 3.625 /im fine-structure line identified by 



Dinerstein fc Geballd (120011 ). This line has recent ly been observed in ab out a dozen PNe be- 



longing to the Milky Way and its dwarf satellites (jPinerstein et al.ll2007l . Dinerstein, Geballe, 
& Sterling, in preparation, hereafter DGS08). Two of these objects exhibit Type I com posi- 
tions: LMC SMP 62, discussed above; and the Galactic PN M 1-40 JGornv et al.ll2004l ). The 
gaseous Zn abundance in M 1-40 is ~ 80% solar, which is consistent with its approximately 
solar Ne and S abundances. Even allowing for modest depletion into dust, which would raise 
the total Zn abundance by 0.1-0.2 dex (assuming conditions typical of the warm ISM; see 
§1), it does not appear that Zn is significantly enriched in M 1-40. In LMC SMP 62, the 
measured Zn ab undance is lower relativ e to solar than the abundances of a elements such 
as O, Ne, and S (IDinerstein et al.l 120071 ) by a factor very similar to the deficiency in Se. In 
this context, it is interesting that several of the iron-group elements, including Ni and Cu 
(which are adjacent to Zn), display behavi or similar to light n-capture elements, which are 
also deficient relative to Fe in LMC giants (jPompeia et al.ll2006l : IMucciarelli et al.ll2008l ). 



3. The stellar models 



The numerical metho d we use has previously been described in detail (ILugaro et al. 



2004J : iKarakas et al.l 120061). We first compute the stellar structure using the Mt Stromlo 
Stellar Structure code (jLattanzio 1986 ). where each model is evolved from the zero-age main 
sequence to near the tip of the TP-AGB. For the models considered in this study, we do not 



^Using Spitzer spectra, Bernard-Salas et al. ( 2008h find higher Ne/H than Leisv fc Dennefeld ( 2006 ). 
However, they adopt a higher value for the solar Ne abundance, and consequently their ratio to solar is 
similar. They also derive a low value for S/H, but the significance of this result is unclear, since these 
authors also find that Galactic PNe have subsolar sulfur abundances ( Pottasch fc Bernard-Sala^ 20061 ). 
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evolve to the post-AGB phase due to convergence difficulties that o ccur when the envelope 



mass is reduced below ~ I.SMq; see iKarakas fc Lattanzid (120071 ) for more details. The 



structure is then used as input into a post-processing nucleosynthesis code where we obtain 
abundances for many more species (up to 166) than are included in the structure model (6 
species). 

The stellar models we computed (see Table [1]) cover a mass range of 3.O-6.5M0 and 
metallicities 0.2-1.0 Zq. The metallicities were chosen to reflect the composition of most 
of the PNe in the SD08 sample, where O and Ar abundances were used as metallicity 
indicators. We include a 2.5Mq, Z = 0.008 model as an example of a low-mass, low- 
metallicity progenitor t hat would lead to a non-Typ e I PN composition. All models use scaled 
solar abundances from lAnders &: Grevessd (Il989l ) except for the Z = 0.012 models, which 



are co mputed with the r evised solar elem ental abundances from lAsplund. Grevesse. &: Sauval 



(120051 ) for comparison. iReimerd ( 

the parameter r] = 0.4, and the 
comparison, we include results for a 5Mq, Z 
loss on the AGB, with the parameter r] = 3.5. 



975r) mass loss was used o n the first giant branch, with 



Vassiliadis fc Wood! (119931 ) mass loss on the AGB. For 



0.02 model computed with Reimers mass 



In Table [T] we summarize st ructural information about t he models. Most of the struc- 
ture calculations are ta ken from iKarakas fc Lattanzid (120071 ). with the Z = 0.012 models 
discussed previously in IKarakas et al.l (120071 ). In Table [1] we present the initial mass and 
metallicity, Z; the total number of TPs computed; the maximum temperature in the He- 
shell, T^^^] the maximum temperature at the base of the convective envelope, Tbc^^; the total 
mass dredged into the envelope during the TP-AGB, Mdred; the final envelope mass Mem', 
and whether or not HBB occurs. In Table (H masses are given in solar units, temperatures 
in millions of kelvins, metallicity Z is the fraction by mass of metals, and elemental and 
isotopic ratios are in terms of number ratios. 



3.1. Type I planetary nebulae and HBB 



Observations suggest that the Typ e I PN progenitor mass range is ~ 3 to 8 Mq (e.g., 
Peimbertlll990l : IStanghellini et al.ll2006l ). For the purpose of this study we consider models 
with masses between 3 to 6.5 Mq, noting that more massive AGB stars would likely evolve 
too quickly to form observable PNe. This mass range encompasses the full range of AGB 
nucleosynthetic behavior: (1) models that experience the TDU but not HBB, and (2) models 
that undergo the second dredge-up (SDU) and HBB in addition to the TDU. From inspection 
of Table [H, the 2.5, 3, and 4 Mq models fall into the first category while the 5, 6, and 6.5 Mq 
models fall into the latter category. The more massive objects that experience both the SDU 
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and HBB are often defined as intermediate-mass AGB stars (see lKarakas et al.ll2006l : iHerwig 
2005h . 



Only models with HBB (i.e., models with initial masses m > SMq) would evolve into 
Type I PNe. Hence efficient extra mixing between the base of the convective envelope and 
a region hot enough to allow some CN cycling is presumably operating to produce the 
observed Type I nucleosynthethic signature of high He/H and N/0 ratios in the lower- mass 
~ 3-4Mq stars. The physical mechanism that could cause such efficient extra mixing is not 
known, although stellar rotation is one possible candidate (see below). It i s probably not 



thermohaline mixing which is on ly efficient in stars less massive than 2Mq (lEggleton et al. 



20061 : ICharbonnel &: Zahru 120071 ). The idea that extra- mixing processes are operating in 
low-mass carbon-rich AGB stars is supported by observation al evidence, includ ing lower 
than predicted ^^C/^'^C ratios observed in N-type AGB stars (lAbia fc IsernI 119971 ). The C, 



O, and Al isotopic ratios found in meteoritic silicon carbide (SiC) and oxide grains from 



(Busso et al. 


1999; 


NoUett et al. 


2003; 


Z inner et al. 


2006; 


Z inner 


2008) 



mixing should not be so efficient to prevent the formation of a carbon-rich atmosphere. 

Some fractio n of Type I and bipolar PNe presumably formed from low -mass stars via 
binary evolution JSokerlll997l : iBahck fc FranMbooi iMoe fc De Marcoll2006l ). although it is 
unclear how this would lead to the high He/H and N/0 ratios observed. Two-dimensional 
simulations of a star with a point-mass companion in a very close orbit suggested that close 
binary evolution has little ef fect on the deep interior s tructure, despite producing large dis- 
tortions to the outer layers (jPeupree fc Karakasll2005l ). However, these simulations did not 
go beyond the point where the stars begin to fill their Roche lobes. Rapid stellar rotation in 
combination with magnetic fields has also been proposed as a mechanism to shape bipolar 
PNe and to produc e the overabundances of He and N in single intermediate-mass progen- 



itors (M > 3Mr^) fICalvet fc PeimbertI Il983l ; iGornv et all 119971 ; iGarcfa-Segura et all Il999 



Dobrincic et al. 20081 ). Further modeling is required to address how effectively binary inter- 
actions, rotation, and/or magnetic fields can affect the internal structure and nucleosynthesis 
of AGB stars. 

In Fig. [3] we plot the temporal evolution of the core, and the temperature in the He- 
shell and at bottom of the convective envelope for the 6.5Mq, Z = 0.012 model. The mass 
of the H-exhausted core shows evidence of efficient TDU whereas the He-shell temperature 
shows that the majority of the TPs reach peak temperatures greater than 300 xlO^K, the 
temperature at which the ^^Ne neutron source is activated. HBB prevented several models 
in Tabled] from becoming carbon rich, and also resulted in low ^^C/^^C ratios that approach 
the equilibrium ratio of ~ 4. The He/H and N/0 ratios are also higher for models with 
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HBB compared to the lower mass models. HBB is more efficient at lower metallicities for a 
given mass owing to the higher temperatures reached at the base of the convective envelope. 
For example, the 6.5 Mq, Z = 0.012 model (corresponding to the Asplund et al. solar 
abundances) displays a higher N/0 ratio than the Z = 0.02 model of the same mass. The 
final N abundance is similar in each of these two models, despite a difference of a factor 
of 1.86 in the initial N abundances, showing that CNO processing is more efficient in the 
low-met alii city case. 



3.2. The nuclear network 



The most important addition to our models for this study is the extension of the nu- 
clear network to include elements heavier than iron up to niobium (M. van Raai et al., in 
preparation). We h ave also made an upd ate of our reaction-rate library - starting from the 
library d escribed in iKarakas et al.l (120071 ) we have included neutron capture cross sections 
from the Bao et al. (j2000 ) compilation. We include 166 species from protons to sulfur and 
iron through to Nb, with a total of 1285 reaction rates corresponding to all /3-decay, p, a, 
and n-capture reactions on all species in the network. We also have two other networks 
that we use to obtain information for the S and Ar isotopes. The ffist has 125 species and 
includes all stable species from H to ^^Ni, while the s econd has 156 species, and includes all 
stable species from H to arsenic (IKarakas et al.ll2007l ). 



As done in lKarakas et al.l (120071 ). we use a "double neutron-sink" description (IHerwig. Langer. fc Lugar 



20031 ) to account for neutron captures on species heavier than that at the end of the network 
(^^Nb in this case). The two artificial species are linked by the reactions ^^Nb(n,7)^^g and 
^^g(n, L)^^g, where ^^g has an initial abundance equal to the sum of the solar abundances 
from Mo to Bi. The second artificial particle, L, is equivalent to the number of neutrons 
captured beyond Nb. The ratio (L/^^g) is a description of the neutrons captured per seed 
nucleus and could in principle be related to the s-process distribution. Note that in the 156 
and 166 networks, the species ^^S is the sum of all species from ^^S to Mn, and the reaction 
^^S(n,7)^^S is assigned an averaged cross section value in order to represent all nuclei from 
3^8 to Mn. 



3.3. Neutron sources in intermediate-mass AGB stars 



The s process is driven by the production of free neutrons that are subsequently cap- 
tured by Fe-peak seed nuclei to form heavier elements. The ^^Ne(a, n)^^Mg reaction was 
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first suKKested as a neutr on source in stars by lCameronI (Il960l ). Later. iTruran fc IbenI (119771 ) 
and ICosner et al.l (119801 ) suggested it was the dominant source in intermediate-mass AGB 
stars, because the He-shells of these stars reach high enough temperatures (T > 300 x 10^ K) 
to allow for this reaction to be efficiently activated (see the middle panel of Fig. [3] as an 
example). The s e tem peratures are reached only in the last few TPs of a lower mass star. 



Truran fc IbenI (119771 1 found that the ^^Ne source results in enhanced levels of the n-rich ele- 
ment Rb owing to the high neutron densities (up to ~ 10^'^ neutrons/cm^), as well as increases 
in the heavy Mg isotopes ^^Mg and ^^Mg from the competi n g reactions ^^Ne(a, n )^^Mg and 



22Ne(a,7)26Mg (iKaeppeler et al.lll994l : iKarakas et al.l 120061 ). tenner et all ( l2003h compared 
results from a chemical evolution model with observations of the neutron-rich Mg isotopes 
and concluded that an extra production site besides Type II supernovae was necessary, and 
that massive AGB stars are a good candidate for such a site. There is, however, some uncer- 
tainty at what Galactic epoch AGB stars started contributing to the chemical enrichment 
of the Galaxy (jSimmerer et al.l l2004l : iMelendez fc CohenI 120071 ) . It is also unknown to what 
exte nt intermediate-mass stars contribute to the Galactic inventory of s-process elements 



[e.g. 



Travaglio et al.ll2004[ l. 



The other potential source of neutrons in AGB stars is the ^^C(a, n)^^0 reaction, which 
operates at lower temperatures (T > 90 x 10^) than the ^^Ne source. Observational and 
theoretical evidence suggests this is th e dominant neutron source in low-mass AGB stars 
( ISmith et al.l 119871 : iGallino et al.lll998l ). To operate efficiently, this reaction requires more 
^^C than is left over from CN cycling in the H-shell; hence some mechanism to mix protons 
from the H-rich envelope into the intershell is needed to produce the extra ^^C. In our 
models, protons are mixed into the intershell region by artificially adding a partial mixing 
zone (PMZ) at the deepest extent of each TDU. These protons are quickly captured by 
the abundant ^^C to form ^^C and ^^N, resulting in the formation of a ^'^C pocket. In the 
^^C pocket, neutrons a re hberated by the reaction ^^C(a, n)^^0 during the interpulse period 
( IStraniero et al.lll995l ). in contrast to the ^^Ne source, which operates during TPs. The 
timescales for neutron production during the interpulse are much longer (> 10^ years) than 
during the convective pulse (~ 10 years), resulting in much lower neutron densities (~ 10^ 
neutrons/cm^) than the ^^Ne source. Together, the timescale for neutron production and 
the neutron source determine the resulting s-process element distribution. The details of 
how the ^^C pocket forms and its extent in mass in the He-intershell are still unknown, 
although various mec hanisms have been proposed, including convective overshoot, rotation, 
and gravity waves; see lHerwigi (120051 ) for a discussion of the relative merits of each mechanism. 

The importance of the ^^C pocket and the ^'^C(q;, n)^^0 reaction for nucleosynthesis 
in intermediate-mass stars is neither clear nor well studied. Adding a PMZ artificially 
into massive AGB models has little effect on the nucleosynthesis of nuclei lighter than Fe 
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(iKarakas et al.ll2006l ) and on the production of Ge in 5Mq, Z = 0.02 models (iKarakas et al. 
20071 ). The first point can be understood by noting that the ^^C pocket is about ~ 10% of 
the mass of the He-intershelL Since the He-intershell is aheady smaller in massive AGB stars 
by about one order of magnitude compared to lower-mass stars, the inclusion of the PMZ 
has little effect on the surface abundances of light elements, which are relatively abundant 
in comparison to heavy elements produced via the s process. The effect of the ^^C pocket 
on the production of s nuclei is less clear, given their low initial abundances. Se and Kr 
elemental abundances derived from Type I PN spectra, in combination with observations 
of the elements Rb and Zr in massive AGB stars in the Galaxy, LMC, and SMC, may pro- 
vide a way of distinguishing between the relative importance of the two neutron sources in 
intermediate-mass AGB stars. 

For the 3, 4, 5, and 6.5 Mq solar composition models we present results with and without 
a PMZ (see Table [2] for the PMZ masses used in each model). No PMZ was included in any 
of the low-metallicity models, with the exception of the 2.5 Mq, Z = 0.008 model, where we 
only include results with a PMZ. The PMZ masses were chosen such that the resulting ^^C 
pockets were ~ 10% to 15% of themass of the He-intershell, consistent with the ^^C pocket 
m asses used by iGallino et al.l (119981 ). The method we use to include a PMZ is described 



m 



Lugaro et al.l (120041 ) , and is similar to that used by iGoriely fc Mowlavil (120001 ) . At the 



deepest extent of each TDU episode, we add an exponentially decaying proto n profile that 



covers the mass specified in Table [2l This method is different from that used by lGallino et al. 
(119981 ). who include a ^^C profile directly into the intershell. 



4. Model Results 

4.1. Germanium to Strontium 

Surface abundances at the tip of the TP-AGB for Zn to Sr are presented in Table ^ 
and Figure H] depicts enrichments of Se and Kr as a function of pulse number in selected 
models. At solar metallicity, the 3Mq models with a ^^C pocket produce the largest surface 
enrichments of s-process nuclei, with final [Se/Fe] and [Sr/Fe] abundances of ~ 0.50 and 0.85 
for the Z = 0.02 model, and 0.56 and 1.2 for the Z = 0.012 model. In the lower metallicity 
models, the largest surface enrichments are seen in the 2.5M0, Z = 0.008 model with final 
[Se/Fe] and [Sr/Fe] abundances of 0.60 and 1.3, and in the 5Mq, Z = 0.004 model with 
~ 0.90 and 0.93, respectively. 



The solar metallicity intermediate-mass (M > 5Mp)) mode l s exh ibit little s-process 
enrichment ([Kr/Fe] and [Sr/Fe] < 0.25). Using lAsplund et al.l (120051 ) solar abundances. 
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the 6.5M0, Z = 0.012 model shows shghtly larger enhancements ([Kr/Fe] = 0.40, [Sr/Fe] 
= 0.30), albeit only when a PMZ is added. The inclusio n of a ^^C pocket ha s little effect 



(< 0.1 dex) on the production of Se and Ge, as noted by iKarakas et al.l (120071 ). but it does 
enhance the Kr and Sr abundances. Only the lower-metallicity 5 and QMq models produce 
substantial enrichments of s-process elements. These models experience a large number of 
TPs (> 60), each with efficient TDU, resulting in significant surface enrichments. No ^^C 
pocket was used for these models, and thus it is clear that the higher temperatures obtained 
during TPs efficiently activates the ^^Ne neutron source. Moreover, the number of neutrons 
per seed nucleus is larger at lower metallicity owing to the smaller abundance of seed nuclei 
(e.g., ^^Fe). 

We find that Kr is enhanced in the s tellar envelope s to a greater degree than Se, in 



agreement with the observational results of ISterling et al.l (120071 ) and SD08. This is because 



Kr lies nearer to the first s-process peak than Se, and has an isotope (^^Kr) with a magic 
number of neutron^. Furthermore, the larger s-process enhancement of K r relative to Se 



i n the AGB models are in agreement with the chemical evolution results of iTravaglio et al. 



( 120041 ) ■ who find that the main s-process contributions of AGB stars to the Solar System Se 
and Kr abundances are 14% and 30%, respectively (the weak s-process, operative in massive 
stars, contributes 25% and 20% to the solar Se and Kr abundances, while the r-process 
provides the remainder). 



The [Ge/Fe] abundances in Tableware different from those given in Table 2 in lKarakas et al. 



( 120071 ). owing to small changes to the neutron-capture cross sections used in the 166 and 
156 species networks. There is a 0.08 dex change in the [Ge/Fe] abundance for the 3Mq, 
Z = 0.02 model with PMZ of O.OO2M0. This difference of 0.08 dex is less than found when 
considering other model uncertainties. For example, there is a 0.16 dex difference in the 
[Ge/Fe] abundances from the two different 5Mq, Z = 0.02 models with different mass-loss 
laws. 

Efficient activation of the ^^Ne neutron source leads to positive (5^^Kr/^^Kr ratio^, as 
can be seen in Table |2] for the most massive AGB models. These positive 5'^^Kr/^^Kr ratios 
are a result of neutron densities in excess of 10*^ cm~'^, which enables neutron captures on ^^Kr 
that produce ®^Kr. If the density is lower, as is the case when ^^C(a,n)^^0 is the neutron 



^Nuclei with a magic number of neutrons are relatively stable against n capture due to their small n- 
capture cross sections, and are the cause of the peaks in the s-process enrichment distribution. 



The S notation is commonly used in meteoritics (see, e.g., iLugaro et al.ll2003af ). and is a variation from 
solar in parts permil, or parts per thousand (where a 100 permil change in the isotopic value is equivalent 
to a 10% variation). 
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source, the branching is bypassed and ^^Rb is produced instead, resulting in negative values 
of the 5^^Kr/*^Kr ratio. The Z = 0.02 models reflect this trend, showing an increase in the 
5'^^Kr/^^Kr ratio with mass, with SMq the minimum mass for efficient activation of the ^^Ne 
source. The ratio is negative in the 3-4Mq models, where the ^^C reaction dominates as the 
neutron source and the ^^Ne reaction is only marginally activated. Values of 5(^^Kr/^^Kr) 
measured in meteoritic stellar Si C grains frorn AGB stars range roughly from —500 to 900, 
and increase with the grain size (ILewis et al.lll994l ). Negative values are well explained by 
AGB model predictions, however, it is difficult to explain the large observed positive values 
because they are only produced by HBB models, where the C>0 condition necessary for the 
formation of SiC grains is not attained, except in the lowest metallicity, Z = 0.004 models. 
Further detailed work is required on this issue. 

The resulting elemental abundance pattern also differs, as Sr is produced in favor of Rb 
in the low-neutron density case, resulting i n low Rb/Sr ratios (~ 0.05). This is c onsistent 
with observations of low-mass AGB stars (ILambert et al.lll995t lAbia et al.l l200ll ). In the 
high-neutron density case, Rb is produced instead, as observed in massive 0-rich AGB stars 
( iGarcia- Hernandez et al.l 120061). Prelimina r y resu lts for Zr and Rb from intermediate- mass 



van Raai et al.l (120081 ). and will be presented in more detail in a 



AGB stars are discussed in 
forth-coming paper (M. A. van Raai et al., in preparation) 



4.2. The iron-peak elements Iron and Zinc 

In Table m we show the [Zn/Fe] ratios and AFe values for each model, where AFe is the 
logarithmic change to the surface iron abundance from the main sequence to the tip of the 
TP- AGB. The AFe values are all very nearly zero, indicating that the surface abundance 
of this element is unaltered by AGB nucleosynthesis. The maximum effect among the Z = 
0.02 models is 2%, and even for the Z = 0.004 case the Fe abundance decreases by less than 
5%. This makes Fe a suitable reference element for estimating metallicities in AGB stars. 

Although the surface Fe abundance is not significantly changed by nucleosynthesis dur- 
ing the AGB, the abundance of Fe in the He-intershell can be depleted by up to ~ 1 dex in a 
massive AGB star (e.g., 6.5 Mq model). In the intershell ^^Fe is destroyed by neutron cap- 
tures, and as the dominant isotope it primarily determines the elemental abundance (even 
if the abundances of the rare neutron- rich isotopes ^^Fe and ^^Fe increase). In lower mass 
models (e.g., 3 Mq), depletion of Fe is restricted to the ^^C pocket that encompasses 10 to 
15% of the mass of the intershell. The next TP homogenizes this small, Fe-depleted region 
with the muc h larger Fe abundance in the rest of the shell. Our 3 Mq results are consistent 



with those of IWerner fc Herwig (120061 ) . who found that AGB models cannot account for 
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the large Fe depletions fou nd in some central stars of PNe, parti c ularly the hot, H-d e ficient 
stars of the PG1159 class jMiksa et alibooi iRauch et aljboosh . IWerner fc Herwiej J2006h 
estimated that these central stars have mean masses of ~ O.62M0, and thus represent the 
evolutionary remnants of lower-mass progenitors than those addressed here, and should be 
compared with models of stars with < SM©. Note that similar or even more extreme Fe 
deficiencies - along with dramatic overabundances of n-capture elements - are seen in older 
white dwarfs where mechanisms such as gravitational settling and levitation b y selective ra- 



diation pressure play irnporta nt roles in establishing the surface abundances (IChayer et al. 
2OO5I : IWerner fc Herwiej Eooef ) . 



In PNe, Fe can be depleted into dust by large factors that may vary from o bject to object 



( Perinotto et al.lll999l ). and even in different regions within the same nebula (ISterling et al. 



20051 ). so that its gas-phase abundance is not a reliable indicator of the true elemental 
abundance. The heavier iron-group eler nent Zn exhibits m uch milder depletions in the ISM, 
particularly in the warm phases (e.g., IWelty et al.lll999l ). Consequently, Zn is frequently 
used as a proxy for iron and the iron-group, for example in damped Lyman-a absorbers 



(e.g.. lAkerman et al.ll2005l ). However, the use of Zn as an indicator of the initial composition 
of an AGB star or PN requires that the change in the abundance of Zn due to internal 
nucleosynthesis is small compared to its initial abundance in the star. It can be seen from 
Table [2] that the predicted changes in the Zn abunda nce are modest o r negli gible for most of 
the models. This is consistent with the conclusions of iTravaglio et al.l (120041 ) . w ho estimated 
the to tal contribution of AGB stars to the solar Zn to be ~ 3%, and with ICayrel et al. 



(l2004j )'s finding that Zn production in the Milky Way is dominated by other sources (i.e., 
explosive Si burning in massive stars). However, at low metallicities, the models show a large 
enough increase that the Zn enrichment may be detectable. 

These results can be compared to the Zn abundances derived for two Type I PNe 
(DGS08). Of the models presented here, the Galactic disk PN M 1-40 is most closely matched 
by the 6.0 M(Tt mod el with Z = 0.02, given its solar Ne and S abundances and N/0 = 1.13 
(iGorny et al.ll2004l ). This model predicts a very small enrichment of Zn, [Zn/Fe] ~ 0.038, 
which is consistent with the measured value (see § 2), even allowing for a possible small 
depletion effect. The case of LMC SMP 62 is more problematical. According to its a-element 
abundances, the most appropriate models are those for Z = 0.008. At this metallicity, all 
three of the models in Table [2], calculated for 2.5, 5.0, and 6.0 Mq respectively, predict 
similar Zn enhancements of ~0.15-0.2 dex. However, as described in § 2, the measured Zn 
abundance in this PN is lower, rather than higher, than those of the a-elements (e.g., Ne, S), 
in contradiction to the models. The disagreement is even worse for Se, which is predicted to 
be enhanced by 0.4-0.6 dex, depending on the initial mass. A modest enrichment in both Zn 
and Se is possible, if their initial abundances were low enough relative to Fe, but nevertheless 
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it remains difficult to satisfy the observational constraints for LMC SMP 62 with the listed 
models. We note, however, that these particular models also predict N/0 values either lower 
(the 2.5 Mq model) or higher (the 5 and 6 Mq models) than the observed value for this 
PN; in § 5 we discuss other possible scenarios for explaining the observed properties of LMC 
SMP 62. 



4.3. Neon through to Argon 



Abundances of S, CI, and Ar in PNe are powerful tools for studying the chemical 
evolution of galaxies and stel lar systems, under the assumptions that these elements are 
unaltered by AGB evo lution (IDopita et al.l 119971 : IStasihska et al.l 1 19981 : iLeisy fc Dennefeld 



2006 



Idiart et a. 



( IStasihska et al 



20071) and are produced in the same relative proportions in massive stars 
19981 ). The case of Ne is not as clear, although the remarkably constant 



Ne/ O ratio of ~ 0.20 (by mass) observed in PNe with a wide variety of chemical compositions 
indicates that in most cases this ratio is not altered during the AGB. ^°Ne, the most abun- 
dant isotope of Ne, is not altered significantly by H or He-burning in AGB stars, whereas 
^^Ne can be produced by successive a-captures onto ^^N. Therefore an enhancement in the 
Ne/0 ratio is expected only w hen ^^Ne is sufficiently enriche d that its abundance is compara- 
ble to or exceeds that of ^°Ne. iKarakas fc Lattanzid (l2003bl ) predicted that such an increase 
would only occur in a narrow range of progenitor masses, ~ 2.5 to ~ 3.5Mq depending on 
metallicity. Our current models are generally consistent with these earlier results, although 
from Table [3] it is clear that the 6.5 Mq, Z = 0.012 model and the 5 and 6 Mq, Z = 0.008 
and 0.004 models have enhanced Ne/0 ratios as a result of some O destruction by HBB, and 
mild enrichments of Ne via dredge-up (maximum Ne increase of ~ 0.3 dex for the 5, 6Mq 
metal-poor models). Therefore, Type I PNe are not expected to be significantly enriched in 
Ne except at low metallicities. 

The ^^Ne is produced primarily in the TPs themselves, but in lower mass models (2.5 and 
3Mq), the partial mixing of protons can further enhance the Ne/0 ratio. This is because the 
partially mixed zone results in a pocket of ^^C and ^^N, thus increasing the ^"^N abundance 
in the intershell to values higher than that from CN-cycling during the interpulse. The 
next TP quickly converts the ^^N into ^^Ne. The resulting ^^Ne abundances after the TP 
are slightly higher than in the case with no partially-mixed zone. After many TPs and 
TDU, the PMZ resu l ts in an increase in the elemental Ne abundance. The observations by 
Bernard-Salas et al.l (120081 ) are intriguing in this respect. They determined Ne/S ratios in a 
sample of 25 PNe in the Magellanic Clouds, and found that four of the objects show high 
Ne/S ratios. For two of the objects, they attribute the high ratios to low S abundances, but 
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in the other two PNe, the high Ne/S values could be due to Ne enrichments. 



Table E] shows the surface abundance predictions for S, CI, and Ar. These elements 
generally are not altered by AGB nucleosynthesis processes by more than the observationa l 
errors (~ ±0.1 — 0.15 dex), consistent with the observations of PNe by lStasihska et al.l (119981 ). 
Neutron captures in the He-shell can result in variations in the isotopic ratios involving the 
stable CI and Ar isotopes; this is demonstrated in the last column of Table [3] where we 
show the permil variations for ^^Ar/^^Ar. Values of 5(^^Ar/^^Ar) meas ured in stellar SiC 
grains range from zero (solar composition) to positive values up to ~ 2 00 (ILewis et al.lll994j ). 
in qualitative agreement with our model predictions. iKahane et al. compared AGB 

model predictions to measurements of CI, Mg, Si, and S isotopic ratios in the cool carbon 
star IRC+10216, and inferred an initial mass of < 2Mq based on the accuracy of the CI and 
Mg data. Kahane et al. concluded that a larger mass progenitor would lead to too high CI 
and Mg isotopic ratios at the surface, although our models suggest that progenitor masses 
as high as ~ 3Mq produce a good fit to the observed data, within the errors quoted for the 
isotopic ratios. In the case of our low-metallicity massive AGB models, neutron captures 
result in modest increases (~ 0.25 dex) in the elemental CI surface abundance. 

The phosphorous abundance can be determined in hot post-AGB objects such as the 
PC 1159 stars, whose envelopes are composed of He-intershell material. Observations in- 
dicate that PC 1159 stars have solar P abundances, whereas 3Mq models predict P abun- 
dances ranging from four times solar to ~ 25 times solar in the He-intershell, depend- 
i ng on the amount of ex tra mixing included in the calculations to produce the ^^C pocket 
( IWerner fc Herwigll2006l ). Our He-shell predictions are consistent with Herwig's model, in 
that we estimate P enhancements from 1.3 to ~ 9 times solar for the 3 Mq model, without 
and with a PMZ respectively. The large P abundances in the He-intershell abundances do 
not result in large surface enrichments (^0.1 dex in all solar metallicity models), so that 
if measurements could be made of the P abundances in PNe, we would not expect to see 
variations from solar. 



Sulf ur is also observed in PG 1159 stars, with abundances varying from 0.01 times solar 
to solar (jWerner &: Herwigl l2006l ) . It is not clear if this range reflects the initial abundance 
of the star, various degrees of nuclear processing, condensation of S into dust or problems in 



determining accurate S abundances. Note that S abundances in several Galactic (IHenry et al. 



20041 ) and some Magellanic Cloud PNe (IBernard-Salas et al.l 120081 ) are low relative to O. 
While infrared observations of S^^ in Galactic PNe seem to indicate that uncertainties in 
corrections for unobserved ionization stages are not the reason for the low S abundances, 
it is possible that S may be d epleted into dust in some PNe (especially C-rich objects; 
Pottasch fc Bernard-Salad I2OO6I ) . AGB models do not predict large depletions of S in the 
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He-shell, with S abundances varying between 0.6 to 0.9 times solar (jWerner fc Herwig||2006l ). 
Our values lie from no change to ~ 0.9 times solar. The results for P and S are puzzling, 
indicating on the one hand that a smaller PMZ is required (to produce less P), or, on the 
other hand that more mixing is required (to destroy more S). 



4.4. Oxygen in intermediate-mass AGB stars 



Oxygen constitutes a large fraction of the total metal content of a star, and in this 
sense is a more fundamental measure of metallicity than Fe/H, which historically has been 
used as the primary indicator of metallicity due to its ease of measurement . While it can 
be challenging to determine O abundances in stars (e.g.. iRamfrez et al.ll2007l ). 0/H is easily 
measured in ionized nebulae because of its bright emission lines, and the resulting abun- 
dances are considered to be among the most accurate for any element. According to the 
current paradigm for galactic chemical evolution, O production is dominated by massive 
stars (M > 10 Mq). Nucleosynthesis models of AGB stars i ndicate that, collectively, these 
stars do not contribute sig nificantly to the production of O fiForestini &: Charbonnell 119971 : 
Karakas &: Lattanzidl2007l ) , although there is evidence for dredge-up of self- produced oxygen 



in a few special cases, includ ing PNe from metal-poor progenitor stars (e.g.. |Pequignot et al. 



2000l : binerstein et allboosh . The inclusion of diffusive convective overshoot at the mner 
edge of the flas h-driven conv ective zone during a TP can lead to O abundances in the inter- 
shell of ~ 20% (jHerwie 2000 ) , i n contrast to standa rd models that predict 2% by mass or less 
( iBoothroyd fc Sackmanrull988l : iKarakas et al.ll2002l ) . This can lead to increases in the surface 
composition of O, particularly at low metallicities. Evidence for this ove rshoot comes from 
PG 1 159 stars that show compositions consistent with Herwig's models ( IWerner fc Herwig 
20061 ). It is unclear however, if this overshooting occurs for all AGB stars or only for AGB 
stars that produce PG 1159 stars as a consequence of late and very late TPs. For exam- 
ple, the observed abundance analyses in intrinsic an d extrinsic AGB sta. rs, as well as in SiC 
grains, suggest that such an overshoot is uncommon (ILugaro et al.ll2003bl ). This is because in 
the region of the He-rich intershell where primary O is synthesized, the temperature is large 
enough for efficient activation of the ^^Ne(Q;, n)^^Mg neutron source, even in low-mass AGB 
stars. In this case, the resulting neutron density is very large (up to ~ 10^^ neutrons/cm^) 
in contrast to observations of low-mass AGB stars that indicates neutron densities of about 



10'' neutrons/cm^ (lAbia et al.ll200ll ). 



Intermediate-mass AGB stars on the other hand, destroy O via HBB when the temper- 
ature at the base of the convective envelope is sufficiently high (T > 80 x 10^ K). The 0-Na 
anti-correlation observed in globular clusters stars has in fact been associated with HBB in 
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intermediate-mass AGB stars, because HBB can qualitatively account for the destruction 
of O along with the production of Na. Moreover, SD08 noted that the Ar/0, S/0, and 
Cl/0 ratios of Type I PNe are approximately a factor of two larger on average than in non- 
Type I objects. Since Ar, S, and CI are not significantly affected by nucleosynthesis in low- or 
intermediate-mass stars (see § 4.3), they interpreted this as evidence for O destruction during 
HBB in Type I PN progenitors. This result had also been found in previous investigations 



pro 

(iMarigo et al.ll2003l : lleisv fc Dennefeldl l2006l : iPottasch fc Bernard-Salad I2006D . 



The amount of O destroyed by HBB depends on the temperature at the base of the 
convective envelope. The temperature, in turn, depends on a number of factors including the 
initial mass and metallicity, with hotter temperatures found in models of increasing mass for 
a given Z, or in models of decreasing metallicity, at a given mass. In addition, the predicted 
efficiency of HBB is sensitive to the treatment of convection and mass loss in the models, 
neither of which are well-constrained empiric ally. In particular, the co nvective efficiency 
affects the luminosity and the mass- loss rate (jVentura fc D'Antonall2005l ). and hence these 
properties are sensitive to the convective model that is utilized. Mass loss determines the 
length of the HBB phase, since HBB does not operate once the envelope mass drops below 
about 1.5 Mq. 

In Table |3] we present the O abundances at the surface of the model star at the tip 
of the AGB. The O abundance remains unaltered in the lower mass models (3 and 4M0, 
Z = 0.02), whereas O destruction takes place in models with HBB. Out of the Z = 0.02 
models, the largest amount of O destruction is seen at G.SMq which has a peak temperature 
at the base of the envelope of 86xlO^K. The effect of metallicity is seen most clearly by 
comparin g the 6.5Mp, models w ith Z = 0.02 and Z = 0.012, with the latter model computed 
assuming lAsplund et al.l (120051 ) C, N, O initial abundances. In the Z = 0.012 model the 
peak temperature is higher, resulting in 0.15 dex destruction of O compared to only 0.09 dex 
in the Z = 0.02 model. Note that the O that is destroyed is all converted into N. The largest 
O destruction of 0.26 dex occurs in the QMq, Z = 0.008 model. Overall, none of the models 
suffered O destruction of more than 0.3 dex. 

It is not clear that the low O abundances relative to S, Ar, and CI in Type I PNe can be 
attributed to HBB. While modeling uncertainties in the treatment of convection and mass- 
loss could in principle allow for more efficient O destruction during HBB than we report, this 
inevitably leads to large increases in the N/0 ratio. This result agrees with the findings of 
Marigo et al.l (120031 ) . whose synthetic models with efficient ON-cycling produced N/0 ratios 
considerably higher than those observed in Type I PNe. Those authors suggested that the 
low O abundances may be explained if the initial metallicities of Type I PN progenitors are 
low. However, this does not explain the approximately solar abundances of S and Ar in these 
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objects. Moreover, our models constrain the metallicities of Type I PN progenitors to be 
Z > OAZq in order to reproduce the observed Se and Kr abundances. This important issue 
remains unresolved. 



5. Discussion 

In Fig. [5] we illustrate the s-process abundance predictions from a selection of the stellar 
models. The gray-hatched box represents the approximate region occupied by Type I PNe 
s-process abundances. Typical observational uncertainties allow the maximum enrichment 
to reach «i 0.6 dex in some cases. The intermediate-mass AGB models are, in general, a 
good match to the observed abundances of Se and Kr in Type I PNe. Models in Table [2] 
that show enhancements of Se and Kr greater than 0.5 dex are the 2.5Mq model, the 3Mq 
models, the QMq, Z = 0.008 model, and the low-metallicity 5Mq, Z = 0.004 model, which 
gives large enhancements for all elements between Zn and Sr. The results for 3Mq with 
a PMZ are only in borderline agreement with the observations, although the results with 
no PMZ are consistent with no s-process enrichments. In either case the final composition 
is carbon-rich, indicating that this model is unlikely to be the precursor of a Type I PN. 
Indeed, the results for the 3Mq models are closer to that found for the 2.5Mq, Z = 0.008 
model which has a distinct, non-Type I PNe composition, with a high final C/0 ratio, and 
considerable overabundances of all s-process elements (> 0.5 dex for Se to Sr). The results 
for the 5Mq, Z = 0.004 model suggests that Type I P Ne in the Galaxy do not descend 



from low-metallicity stars, as suggested, for example, by iMarigo et al.l (120031 ). and that a 
minimum metallicity of ~ 0.4 Zq is required for these PNe. This result is not unexpected, 
given that Type I PNe evolve from more massive stars than other PNe, and hence trace a 
younger population. On the other hand, these model results suggest that Type I PNe in the 
Magellanic Clouds could be more highly enriched in light n-capture elements such as Se and 
Kr than their Galactic counterparts. This is a strong motivation for further observations 
of s-process elements in PNe belonging to the Magellanic Clouds or to other metal-poor 
populations. 

Observations of other trans-iron elements in PNe in addition to Se and Kr could be 
used to constrain the mass and metallicity ranges of Type I PNe. For example, Sr and 
Rb abundances are potentially useful discriminants of progenitor mass. Sr is considerably 
enriched in the 3Mq models ([Sr/Fe] ~ 0.85 — 1.23), but is slightly if at all enhanced in the 
6 and 6.5Mq models ([Sr/Fe] < 0.4 dex). Note that the inclusion of a PMZ in the 6.5Mq, 
Z = 0.012 model increases the level of Sr enrichment expected, from 0.05 to 0.3 dex. 



Of these next two elements heavier than Kr (element number 36), Rb (element 37) 
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is more likely to be measurable in a meaningful way in nebulae than Sr (element 38), for 
the following reasons. First, Sr is refractory and can be strongly depleted into dust grains. 
In fact, based on its condensation temperature (ILoddersi l2003l ). Sr is expected to be in- 
corporated into dust to a greater degree than Fe, whose gaseous abundance is depleted by 
1-2 orders of mag nitude in PNe (jPerinotto et al.l Il999l : [Sterling et al.l 120051 ) and the ISM 
( IWelty et al.lll999l ). Therefore, measurements of the gaseous Sr abundance in PNe are not 
likely to reveal useful information about s-process enrichments. Second, the only emission 
lines of Sr in the optical spectral region originate from relatively highly ionized states that 
are not expected to be abundant, except in a fe w of the highest-excitation PNe. Sr emis- 



sion has been identified in only one PN to date (jPequignot &: Baluteau 



20051 ). On the other hand, based on its condensation temperature (iLodders 



1994 



Zhang fc Liu 



2003D Rb does 



not deplete into dust as readily as Sr, and thus its gaseous abundance should be more repre- 
sentative of the overall elemental abundance. Rb emission has been identified in a number 



of PNe flPequignot fc Baluteaulll994l : IZhang fc Liull2005l : ISharpee et al.ll2007r i and should be 
detectable in PNe with a broader range of ionization than Sr. This discussion suggests that 
Rb is a useful element with which to constrain stellar models to PN observations. 

Furthermore, Rb is observed to be substantially enhan ced in massive Galactic OH/IR 
stars with the maxira um enhancement [Rb/Fe] ~ 2 ± 1 dex (iGarcfa- Hernandez et al.ll2006l ). 
van Raai et al.l (120081 ) have made a preliminary comparison between Zr, Rb, and Li predic- 
tions from 5 to Q.5Mq solar composition models to the observations of the OH/IR stars. The 
best match between stellar models and the observations was obtained with the 6.5Mq model 
computed by 1) not i ncluding a ^^C pocket , and 2) by accounting for extra TPs not modeled 
in detail, as done by iKarakas et al.l (120071 ). The extra TPs and TDU mixing episodes may 
occur because not all the envelope mass was lost when convergence difficulties ended the 
computation. The main uncertainty is the behavior of the TDU efficiency at small envelope 
masse s, with studies genera l ly finding a decrease in the efficiency with decreasing envelope 
mass (jStraniero et al.l 119971 : iKarakas et al.ll2002l ). Further studies on the evolution of the 
TDU efficiency with decreasing envelope mass are needed to help settle this issue. The tab- 
ulated results presented in this study were calculated without the inclusion of these extra 
TPs and can be considered lower limits to the final PN abundances, noting that it is highly 
uncertain how much of an impact the final TPs would have on the the final PN abundances. 

The final few TPs can potentially have a sign i ficant effect on the final abundance of the 
star and its resulting PN. For example, iFrost et al.l (119981 ) noted that intermediate-mass AGB 
stars may become luminous, optically obscured carbon stars near the end of the TP-AGB, 
when mass loss has removed much of the envelope, extinguishing HBB but allowing dredge- 
up to continue. Under these conditions the envelope mass is signficantly reduced, resulting 
in only minimal dilution of the intershell material that is dredged up. This can allow for 
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large enhancements to the surface composition of C and other He-burning products. In the 
models with HBB included in this study, all experienced a few TPs and TDU episodes after 
the cessation of HBB, with some of the lower metallicity models becoming carbon rich (see 
Tabled], i.e., GMq, Z = 0.008; 5 and 6Mq, Z = 0.004) before convergence difficulties ended 
the computation. Note that it is after this stage that we consider the effect of remaining 
TPs not modeled in detail using a synthetic AGB algorithm. 



In Ivan Raai et al.l (120081 ) the inclusion of these remaining TPs only had a large impact on 



the surface abundances of the 6.5Mq, Z = 0.02 model (out of the 5, 6 and 6.5Mq, Z = 0.02 
models), and only when assuming that efficient TDU occurs at a very small envelope mas^. 
A final [Rb/Fe] ~ 0.9 was obtained in the 6.5Mq, Z = 0.02 model, in agreement with the 
observed Rb abundances of OH/IR stars given the large uncertainties. Along with increases 
in Rb, the [Se, Kr/Fe] abundances were increased to 0.6 and 0.7 dex (respectively), above 
the observed Type I PN limit of 0.3 dex. This model was estimated to have 7 remaining 
TPs, with the last TP occurring with a very small envelope mass of ~ O.OSMq. The small 
envelope mass during the last TDU leads to the situation discussed above: A substantial 
increase in the surface abundances, caused by the mass of the envelope being only ~ 10 
times greater than the mass of the intershell. It is questionable whether such an efficient 
TDU would occur for such a low envelope mass; so if we exclude this last mixing episode, the 
[Se, Kr/Fe] abundances are ~ 0.4 and 0.45, respectively, in agreement with observed Type I 
PN abundances. This analysis may indicate that the OH/IR stars are not the progenitors 
of most Type I PNe and that there are possible observational biases in the sampling of each 
population. That is, the OH/IR population may sample the most massive AGB stars (6- 
7Mq), with the largest [Rb/Fe] (and [Kr/Fe]) abundances, while the Type I PN population 
may come from slightly less massive AGB stars, with lower Kr and Rb abundances. There is 
also the problem of explaining the large Rb abundances in the OH/IR stars, given that they 
seem to require the TDU to operate efficiently at small envelope masses. More observations 
of Rb in OH/IR stars and PNe will help to settle these issues. 

Modeling uncertainties can dramatically affect AGB nucleosynthesis predictions. As 
noted previously, different convective models have a large impact on the nucleosynthesis of 
light elements (e.g., N, O, Na) in AGB models with HBB. Improvements in our understanding 
of convection are desperately needed but this may only come about through mult idimensional 



studies (e.g.. iDearborn et al.ll2006l : iHerwig et al.ll2006l : iMeakin fc ArnettI 120071 ). The mass 



loss law used in the computations is similarly important, because this can determine (for 
a fixed convective model) the number of TPs and mixing episodes. For example, if we 



^Besides the inclusion of the extra TP the models presented here and in Ivan Raai et al.l (j2008l ) are the 
same. 
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artificially cut the number of TPs by half for the 5Mq, Z = 0.004 model (40 TPs instead 
of the computed 80), the final [Kr/Fe] ratio is ~ 0.35 dex instead of 1.17 dex. This simple 
analysis indicates that we need better constraints on the mass loss from low-metallicity, 
massive AGB stars if we are to provide constraints on the metallicity range of Type I PNe. 

Another importan t uncertainty relate s to the unkno wn formation mechanism for ^^C 
pockets in AGB stars (IGallino et al.l Il998l : iHerwigl l2005l ). These pockets are necessary to 
provide neutrons in low-mass stars. It is unclear, however, what role the ^^C neutron source 
plays in more massive AGB stars. We had hoped to ascertain the dominant neutron source 
from comparisons of Type I PN abundances to massive AGB stellar models, but the Se and 
Kr abundances alone are not sufficient. Selenium abundances were not strongly affected by 
the inclusion of a ^^C pocket in the 5 and 6.5Mq models, whereas the abundance of Kr was 
affected but not increased above 0.5 dex. Note that in three Type I PNe of SD08 where 
both Se and Kr were detected, an abundance of [Kr/Se] =0.5 was found, providing some 
support for the inclusion of a ^^C pocket in massive AGB models; this conclusion is, however, 
limited by small number statistics. We predict that PNe should have larger enrichments of 
Kr than Se. While this result is valid regardless of whether or not a ^^C pocket is included, 
larger [Kr/Se] values are predicted in models with a PMZ (see Table 2). Firmer conclusions 
regarding the neutron source operating in intermediate-mass AGB stars may come from large 
data sets of many ri-capture elements, including Br, Rb, and Xe. 

Finally, models with extra mixing may be necessary to reproduce the observed abun- 
dances of some PNe. For example, the composition determined for LMC SMP 62 (see §2 
and §4.2) suggests that it is a Type I PN, but its low Se and Zn abundances are at odds 
with our predicted abundances for intermediate-mass models. A 3Mq, Z = 0.02 model with 
rotation by Cantiello & Langer (private communication) had a N/0 ratio of 0.43 after the 
first dredge up, compared to 0.32 from our model with no rotation, whereas the N/0 value 
of LMC SMP 62 is ~ 0.4 - 0.5. During core H-burning, rotation mixes processed material to 
layers further out than they would otherwise be found. Rapid rotation may also inhibit the 
production of s-process elements during the AGB phase ( iHerwig et al.ll2003l ). consistent with 
the lack of Zn and Se enrichments found in this object. An alternative is that LMC SMP 
62 evolved from a low-mass progenitor of less than 2Mq, which di d not exper i ence efficient 
TDU. Note that the low C/0 ratio of 0.17 found for this object (lAUer et al.l 119871 ) is also 
consistent with a star that did not experience any TDU mixing (suggesting an initial mass 
of < 1.25Mq). It is unclear if rotation can explain the high N/0 ratio in such a low-mass 
star. Other extra- mixing processes (e.g., thermohaline mixing) may need to be invoked in 
this case. 
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6. Conclusions 

Type I PNe have high He/H and N/0 ratios, indicating that they are the descendants of 



i nterm ediate-mass AGB stars with initial masses between ~ 3 to 8Mq. [Sterling &: Dinerstein 



( 120081 ) found that Type 1 PNe exhibit significantly smaller enrichments of Se and Kr (< 
0.3 dex) on average than other PNe. We calculated s-process enrichments in a set of AGB 
models covering a range in mass from 2.5 to 6.5Mq, and metallicity from 0.2Zq to Zq. The 
2.5Mq, Z = 0.008 model was included to show an example of a model that would produce a 
clear non-Type I PN abundance signature, with [Se, Kr/Fe] > 0.5 and C/0 ~ 4. The main 
conclusion is that the results for the 3-6. 5M0 are a good match to the observed abundances. 
The only real exception out of the HBB models is the low-metallicity 5Mq, Z = 0.004 model 
that produced much larger enhancements in Se and Kr than observed. This suggests that 
Galactic Type I PNe do not descend from such low-metallicity objects. 

We also compare calculated abundances for selected intermediate-mass (Ne, P, S, Ar) 
and iron-peak (Fe, Zn) elements to observations of post-AGB stars and PNe. We find that 
the elemental abundances of P, S, CI, Ar, Fe, and Zn are essentially unchanged by AGB 
nucleosynthesis, although there are isotopic shifts caused by neutron captures in the He- 
shell. These results justify use of elements such as S, CI, and Ar in PNe as tracers of 
Galactic chemical evolution. 

It is difficult to reach firm conclusions about the neutron source operating in massive 
AGB stars from Se and Kr abundances in Type I PNe. Certainly it seems that at the least the 
^^Ne source, with efficient TDU, is required to produce enhancements in s-process elements 
in massive Type I PN progenitors. It is less clear whether a ^^C pocket is also required 
in these stars. Increases in the [Kr/Fe] abundance ratio were observed in the models when 
a ^^C pocket was included, but not beyond the amounts observed in Type I PN spectra. 
Obtaining abundances for more n-capture elements, particularly Rb and elements beyond 
the first s-process peak, from a large data set of OH/IR stars and PNe will help distinguish 
among the possibilities. 

Finally, only the models with HBB (M > 5M0, depending on Z) show the high He/H 
and N/0 ratios that define the Type I PN class. Given that 3Mq stars are more common 
than 6.5M(7) stars (according to initial mass function and evolutionary time considerations), 
it may be necessary that another mixing process other than HBB is active in stars of the 
lower-mass range (3-4M0), if these stars do in fact evolve into Type I PNe. Rapid stellar 
rotation in single ~ 3M0 stars may, for example, be able to account for the increased He and 
N abundances in the progenitor stars. The problems of convection and mixing, and of other 
modehng uncertainties (e.g., ri-capture cross sections) requires further study in the context 
of s-process nucleosynthesis in intermediate-mass AGB stars. 
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Table 1: Details of the stellar models. 



Mass 


Z 


TPs 




T-imax 
bee 


M 


Ircd 


Menv 


HBB? 


C/0 


12C/13C 


N/0 


He/H 


3.0 


0.02 


26 


302 


6.75 


8.1( 


-2) 


0.676 


No 


1.27 


108 


0.303 


0.122 


4.0 


0.02 


17 


332 


22.7 


5.6( 


-2) 


0.958 


No 


0.99 


76.5 


0.336 


0.118 


5.0 


0.02 


24 


352 


64.5 


5.0( 


-2) 


1.500 


Yes 


0.77 


7.84 


0.542 


0.133 


5.0 


0.02^ 


37 


368 


57.0 


1.0( 


-1) 


1.922 


Yes 


1.71 


15.5 


0.465 


0.139 


6.0 


0.02 


37 


369 


83.1 


5.8( 


-2) 


1.791 


Yes 


0.38 


10.8 


1.251 


0.150 


6.5 


0.02 


40 


368 


86.5 


4.7( 


-2) 


1.507 


Yes 


0.40 


11.6 


1.205 


0.154 


3.0 


0.012 


21 


307 


7.23 


9.2( 


-2) 


0.805 


No 


2.80 


190 


0.366 


0.117 


6.5 


0.012 


51 


369 


90.0 


6.5( 


-2) 


1.389 


Yes 


0.76 


10.4 


2.503 


0.146 


2.5 


0.008 


27 


302 


5.42 


1.1( 


-1) 


0.664 


No 


4.17 


392 


0.307 


0.109 


5.0 


0.008 


57 


366 


80.8 


1.7( 


-1) 


1.795 


Yes 


0.93 


7.48 


5.937 


0.128 


6.0 


0.008 


68 


374 


89.6 


1.2( 


-1) 


1.197 


Yes 


1.40 


8.80 


5.799 


0.135 


5.0 


0.004 


81 


377 


84.4 


2.2( 


-1) 


1.141 


Yes 


3.62 


10.8 


16.79 


0.138 



"Computed with Reimers mass loss on the AGB. All other models have I Vassiliadis fc WoodI (|l993r i mass loss. 
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Table 2: Surface abundance results for Zn through to Sr, taken at the tip of the AGB from 
the last computed model. 



Mass 


PMZ 


[Zn/Fe] 


[Ge/Fe] 


[Se/Fe] [Br/Fe] 


[Kr/Fe] ^^^j^r/^^j^r^ 


[Sr/Fe] 


AFe"^ 






Z = 


0.02. Anders & Grevesse (1989) C. N and 








3.0 


0.0 


0.001 


0.004 


0.004 < -0.001 


0.006 


-56.05 


0.007 


< 0.001 


3.0 


0.001 


0.071 


0.263 


0.346 0.299 


0.477 


-693.3 


0.700 


-0.00152 


3.0 


0.002 


0.122 


0.397 


0.493 0.433 


0.628 


-774.4 


0.853 


-0.0021 


4.0 


0.0 


0.004 


0.011 


0.009 0.003 


0.010 


-25.82 


0.007 


-0.0016 


4.0 


l(-4) 


0.077 


0.282 


0.365 0.282 


0.519 


-40.31 


0.587 


-0.0026 


5.0 


0.0 


0.012 


0.031 


0.026 0.011 


0.028 


-10.85 


0.011 


-0.0024 


5.0 


l(-4) 


0.024 


0.096 


0.132 0.083 


0.241 


307.7 


0.246 


-0.0024 


5.0= 


0.0 


0.073 


0.188 


0.181 0.107 


0.210 


91.56 


0.085 


-0.0090 


6.0 


0.0 


0.038 


0.106 


0.104 0.051 


0.131 


92.09 


0.055 


-0.0038 


6.5 


0.0 


0.047 


0.133 


0.138 0.064 


0.187 


186.3 


0.083 


-0.0037 






Z -- 


= 0.012. AsDlund et al. r2005) G. N and 








3.0 


0.002 


0.106 


0.434 


0.566 0.473 


0.645 


-747.2 


1.232 


-0.0028 


6.5 


0.0 


0.042 


0.113 


0.111 0.054 


0.136 


55.60 


0.053 


-0.0047 


6.5 


l(-4) 


0.052 


0.164 


0.203 0.122 


0.397 


1171 


0.300 


-0.0048 


Z = 0.008 


2.5 


0.002 


0.148 


0.478 


0.600 0.547 


0.794 


-722.8 


1.31 


-0.0040 


5.0 


0.0 


0.163 


0.380 


0.401 0.302 


0.506 


355.3 


0.264 


-0.0138 


6.0 


0.0 


0.179 


0.432 


0.476 0.324 


0.657 


712.7 


0.388 


-0.0098 


Z = 0.004 


5.0 


0.0 


0.403 


0.813 


0.897 0.734 


1.176 


1178 


0.927 


-0.0216 



"Si^X/^'^^X) notation is defined according to [{'X/'''^X)/{'X/'''^X)q - 1] x 10^ where '°^X is a reference 
isotope. 

''Defined by log(yfinai/yinitiai), where Y is the surface abundance in number fraction. 
^The model has Reimers mass loss on the AGB instead of Vassiliadis &: Wo'odI ( 1993h . 
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Table 3: Abundance results for O, Ne, P, S, CI, and Ar. 



Mass 




Z 


e(0)^ 

V / 


e(Ne) 


Ne/0^ 


e(P) 


e(S) 


e(Cl) 


e(Ar) 


^38Ar/36Ar= 


initial 




0.02 


8.935 


8.103 


0.181 


5.591 


7.277 


5.143 


6.570 


0.000 


3.0 


no 


PMZ 


8.940 


8.279 


0.273 


5.595 


7.299 


5.192 


6.592 


7.735 


4 n 




PMZ 

J. 1VJ.ZJ 


8 933 


8 204 


234 


5 61 3 


7 294 


5 1 90 


6 586 


14 55 


5.0 


no 


PMZ 


8.931 


8.191 


0.227 


5.617 


7.310 


5.208 


6.603 


17.76 


6.0 


no 


PMZ 


8.908 


8.168 


0.227 


5.696 


7.328 


5.237 


6.621 


26.29 


6.5 


no 


PMZ 


8.905 


8.166 


0.228 


5.635 


7.330 


5.212 


6.623 


25.18 


initial 




0.012 


8.659 


7.838 


0.189 


5.358 


7.139 


5.498 


6.178 


0.000 


3.0 


no 


PMZ 


8.656 


8.086 


0.336 


5.393 


7.165 


5.530 


6.205 


17.00 






0.002 


8.666 


8.666 


0.486 


5.522 


7.165 


5.534 


6.204 


41.34 


6.5 


no 


PMZ 


8.563 


7.868 


0.277 


5.410 


7.188 


5.548 


6.227 


61.39 


initial 




0.008 


8.514 


7.675 


0.181 


5.155 


6.849 


4.715 


6.142 


0.000 


2.5 


no 


PMZ 


8.521 


7.846 


0.431 


5.178 


6.870 


4.754 


6.162 


16.37 






0.002 


8.528 


7.985 


0.653 


5.251 


6.867 


4.762 


6.160 


32.81 


5.0 


no 


PMZ 


8.419 


7.817 


0.313 


5.454 


6.899 


4.886 


6.192 


92.67 


6.0 


no 


PMZ 


8.317 


7.781 


0.364 


5.471 


6.908 


4.861 


6.201 


63.44 


initial 




0.004 


8.206 


7.367 


0.181 


4.847 


6.542 


4.407 


5.835 


0.000 


5.0 


no 


PMZ 


8.052 


7.665 


0.512 


5.444 


6.609 


4.633 


5.900 


136.7 



°e(Y) = log 10{Y/H) + 12, where abundances are measured by number, and H is the abundance of hydrogen. 
^The Ne/0 ratio is calculated by mass, to allow comparison to the observations. 
^The S notation was defined previously in Table [H 
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Non-Type I, Non- Bipolar PNe 
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[Se/(0, Ar)] 




-0.5 
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0.5 1.0 
[Kr/(0. Ar)] 



Fig. 1. — Histograms of Se and Kr abundances from SD08, separated into 0.1 dex bins, for 
non- bipolar, non- Type 1 PNe (top); Type 1 PNe (middle); and bipolar PNe (bottom). Data 
are shown only for objects in which Se and/or Kr emission was detected. Type I and (to a 
lesser extent) bipolar PNe exhibit smaller enrichments than PNe with less massive progenitor 
stars. 
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Fig. 2. — The Se and Kr abundances of PNe plotted against He/H (top panels), N/0 (middle 
panels), and C/0 (bottom panels). The triangles correspond to Type I PNe, the boxes are 
non-Type 1 bipolar PNc, and circles represent all other PNe from the sample of SD08. Sc and 
Kr upper limits arc depicted as open symbols. Representative error bars to the abundances 
arc shown in the left hand panels. The observed correlations show that Type 1 PNc generally 
display smaller s-process and C enrichments than other PNe. While many bipolar PNe follow 
the same trends, some are significantly enriched in C, Se, and Kr, suggesting that they arise 
from less massive stars. 
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Fig. 3. — The temporal evolution of the (a) hydrogen (solid line) and helium-exhausted 
(dashed line) cores, (b) the He-shell temperature and (c) the temperature at the base of the 
convective envelope for the Q.5Mq, Z = 0.012 model during the TP-AGB. 
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Fig. 4. — We show the [Se/Fe] ratio (top panel) and the [Kr/Fe] ratio (bottom panel), as a 
function of the thermal pulse number, from a selected number of stellar models. Each point 
represents the surface abundance during the interpulse period. 
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iM^, 'Z = '0.02' 
*) 6.51^^, Z = 0.02 
(□) 5MJ: Z = 0.008 
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Proton number, Z 

Fig. 5. — The surface abundances of Zn, Ge, Se, Br, Kr, and Sr at the tip of the AGB phase 
for a selection of the stellar models. Abundances are shown as [X/Fe] ratios, and are plotted 
as a function of proton number, Z. The 3Mq, Z = 0.02 model includes a PMZ of O.OO2M0 
whereas all other models are without ^^C pockets. The gray-hatched box with a maximum 
at [X/Fe] = 0.3 represents the region of observed Type I PN Se and Kr abundances. Note 
that observational uncertainties allow for [Se, Kr/Fe] up to 0.6. Boxes around the elements 
Se and Kr indicate that for the present study we can only use these elements for comparison 
as large PN abundance data-sets of the other trans-iron elements are not yet available. 



